Family of Hermitian Low-Momentum 
Nucleon Interactions with Phase Shift Equivalence 
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Using a Schmidt orthogonalization transformation, a family of Hermitian low-momentum NN in- 
teractions is derived from the non-Hermitian Lee-Suzuki (LS) low-momentum NN interaction. As 
special cases, our transformation reproduces the Hermitian interactions of Okubo and Andreozzi. 
Aside from their common preservation of the deuteron binding energy, these Hermitian interactions 
are shown to be phase shift equivalent, all preserving the empirical phase shifts up to decimation scale 
A. Employing a solvable matrix model, the Hermitian interactions given by different orthogonaliza- 
tion transformations are studied; the interactions can be very different from each other particularly 
when there is a strong intruder state influence. However, because the parent LS low-momentum 
NN interaction is only slightly non-Hermitian, the Hermitian low momentum nucleon interactions 
given by our transformations, including the Okubo and Andreozzi ones, are all rather similar to 
each other. Shell model matrix elements given by the LS and several Hermitian low momentum 
interactions are compared. 

PACS numbers: 21.60.Cs 



I. INTRODUCTION 

A fundamental problem in nuclear physics has been 
the determination of the effective nucleon-nucleon (NN) 
interaction appropriate for complex nuclei. Typically, 
one starts from a NN interaction, Vnn, constrained by 
the deuteron properties and the empirical low-energy NN 
scattering phase shifts. Several realistic meson models 
0, HH0 f° r Vnn have been obtained in this way, and 
while they all share the same one pion tail, they differ 
significantly in how they treat the shorter distance com- 
ponents. Despite this difference, these models all give 
approximately the same low-energy phase shifts and 
deuteron binding energy. This result clearly manifests 
the main theme of the renormalization group (RG) and 
effective field theory (EFT) approach, namely physics 
in the infrared region is insensitive to the details of the 
short-distance dynamics 0, 0, 0, 0, 0, EO, E3, E3, E3, US - 
ft is thus possible to have infinitely many theories that 
differ substantially at small distances, but still give the 
same low-energy physics, as long as they possess the same 
symmetries and the "correct" long-wavelength structure. 
Since low-energy physics isn't concerned with these high- 
energy details, one should just use an effective theory 
with the short-wavelength modes integrated out. 

Following this RG-EFT idea, a low-momentum NN 
effective interaction, V; oul _fc, was recently developed 
111 111 EE 13 13 While similar in spirit to tra- 
ditional EFT, Vi ow -k is not derived via the usual RG- 
EFT methods; rather, it combines the standard nuclear 
physics approach (SNPA) with EFT making it a more 
effective EFT (MEEFT) [HEEl, as discussed in 0. 
As such, a main step in its derivation is the integrating 



out of the high momentum components of some realistic 
NN potential model Vnn such as those of Ref. 0,HH0]- 
Even though these Vnn models arc quite different from 
each other, it is remarkable that the Vi ow -kS derived 
from them are nearly identical to each other, suggesting 
a nearly unique Vi ow -k 25]. Furthermore, shell model 
calculations using Viow-k have given very encouraging 
results over a wide range of nuclei [HE>[l!l23- Appli- 
cations of Vi ow -k to quasi-particle interaction, superfluid 
gaps and equation of state for neutron matter have also 
been highly successful [l8l ]. 

An important problem in deriving Vi ow -k is how to 
obtain a low-momentum NN interaction which is Her- 
mitian. The Vi n , ,,-k. given by the T-matrix equivalence 
approach flH Il6j| is not Hermitian, and some additional 
transformation is needed to make it Hermitian. There 
are a number of methods used to obtain a Hermitian ef- 
fective interaction, such as those of Okubo j2(|, Suzuki 
and Okamoto |23, and Andreozzi [2^. Which of these 
methods should one use? How different are the Hermitian 
Vi ow -kS given by them? These questions seem to have 
not been investigated. In the present work, we shall study 
the Vi ow -kS given by these methods as well as develop 
a unified method by which a family of phase-shift equiv- 
alent Hermitian low-momentum NN interactions can be 
obtained. It is, of course, important that Viow-k preserve 
phase shifts, and while the non-Hermitian VJ ou ,_fe given 
by the T-matrix equivalence [15l Il6| can be shown to 
preserve phase shifts in a straightforward way, it seems 
to be more involved to prove the phase shift preservation 
by Hermitian Viow—k- Epelbaoum et al. |8j have pointed 
out that phase shifts are preserved for the Okubo low- 
momentum interaction 8], but phase shift preservation 
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for other Hermitian Vi ow — &'s seems to have not been in- 
vestigated - an issue we will also examine. 

To conclude our introduction, we shall briefly review 
the non-Hermitian low-momentum interaction given by 
the T-matrix equivalence approach. Then in Section II 
we present a general method, based on Schmidt orthogo- 
nalization transformation, for generating a family of Her- 
mitian effective interactions, and show that the Hermi- 
tian interactions of Okubo, Suzuki and Okamoto and An- 
dreozzi all belong to this family. In Section III, we study 
our method using a simple solvable matrix model of the 
Hoffmann type 29] ; focusing on the influence of intruder 
states and the difference of the Hermitian interactions 
given by the various methods. In Section IV we present 
a proof that phase shifts are preserved by the Hermitian 
interactions generated by our method. This preservation 
will also be checked by phase shift calculations using dif- 
ferent Hermitian interactions. Finally in Section V we 
present our results for the Vi ow _kS corresponding to the 
various Hermitian interactions, where we show that our 
approach allows us to construct a Vi ow -k which preserves 
the deuteron wave function, in addition to the preserva- 
tion of deuteron binding energy and low-energy phase 
shifts. 

In the following let us first briefly review the Vi ow -k 
given by T-matrix equivalence |l5l Il6| . 
the half-on-shcll T matrix 



One starts from 



T(k', k, k 2 ) = V NN {k', k) + P q 2 dqV NN {k', q) 

Jo 

X T2 5— ^+T(q,k,k 2 ). (1) 

K — q A + iVT 

An effective low-momentum T matrix is then defined by 
T low ^ k (p',p,p 2 ) = Vi ow - k (p,p) (2) 



+P / q 2 dqVi ow - k {p',q) 



1 



—Ti ow _ k (q,P,p 2 ) 



p 2 — q 2 + z0+ 

where the intermediate state momentum is integrated up 
to A. In the above two equations we employ the principal 
value boundary conditions, as indicated by the symbol P 
in front of the integral sign. We require the T matrices 
satisfy the condition 



T{p ,p,p ) = T low - k (p ,p,p ); (p ,p) < A. (3) 

The above equations define the effective low-momentum 
interaction V/ ot „_fc. Using a Q-box folded-diagram 
method [H|33, it has been shown [HE [13 that the 
above equations are satisfied by the solution 
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Q-Q' Q + Q' Q 



Q' Q Q 



where Q-box denotes the irreducible vertex function 
whose intermediate states are all beyond A, and Q' is 
the same vertex function except that it starts with terms 
second order in the interaction. The low-momentum ef- 
fective NN interaction of Eq.(4) can be calculated using 
iteration methods such as the Lee-Suzuki |3fj|. Andreozzi 
[28| or Krenciglowa-Kuo [3l| methods. 

The above Vi ow -k preserves both the deuteron binding 
energy and the half-on-shcll T-matrix of Vjq n (which im- 
plies the preservation of the phase shifts up to Ei a b = 
being the nucleon mass). This Vi ow -k 
is not Hermitian, as indicated by Eq.(4). As we will 
show soon, starting from this Vi ow -k a family of phase- 
shift equivalent Hermitian low-momentum NN interac- 
tions can be obtained. 



II. FORMALISM 

Before presenting our general Hermitization procedure, 
let us first review some basic formulations about the 
model space effective interaction. We start from the 
Schroedinger equation 



{H + V)9 n = E n %, 



(5) 



where Hq is the unperturbed Hamiltonian and V the in- 
teraction. The eigenstates of Hq are 4> n with eigenvalues 
e n . For example, Hq can be the kinetic energy operator 
and V the NN interaction Vjvjv- A model-space projec- 
tion operator P is defined as I I) where d is 
the dimension of the model space. The projection oper- 
ator complement to P is denoted as Q, and as usual, one 
has P 2 = P,Q 2 = Q and PQ = 0. In the present work, 
P represents all the momentum states with momentum 
less than the cut-off scale A. 

A model-space effective interaction V e ff is introduced 
with the requirement that the effective Hamiltonian 
P(Hq + V e /f)P reproduces some of the eigenvalues and 
certain information about the eigenfunctions of the orig- 
inal Hamiltonian (Hq + V). There are a number of ways 
to derive V e ff, but, as indicated by Eq.(4), our effec- 
tive interaction is obtained by the folded diagram method 
[33,l3^| and can be calculated conveniently using the Lee- 
Suzuki-Andreozzi [2^] or Krenciglowa-Kuo [3l| iteration 
methods. We denote this effective interaction as Vls, 
with the corresponding model space Schroedinger equa- 
tion 

P{Hq + V LS )PXm = E mXm , (6) 
where {E m } is a subset of {E n } of Eq.(5) and Xm = 

It is convenient to rewrite the above effective interac- 
tion in terms of the wave operator ui, namely 



(4) 



PV LS P = Pe-"(H + V)e"P - PH P, 



(7) 
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where w possesses the usual properties 

to = QujP; 
Xm = e~ u 9 m ; 

UXrn = Q*m- (8) 

While the eigenvectors ^ n of Eq.(5) are orthogonal to 
each other, it is clear that the eigenvectors \ m of Eq.(6) 
are not so and the effective interaction Vls is not Her- 
mitian. We now make a Z transformation such that 

Z Xm = v m J 

(w m I %') = $rnm>; m,m' = l,d, (9) 

where d is the dimension of the model space. This trans- 
formation reorients the vectors Xm such that they be- 
come orthonormal to each other. We assume that Xm's 
(to = 1, d) are linearly independent so that Z^ 1 exists, 
otherwise the above transformation is not possible. Since 
v m and Z exist entirely within the model space, we can 
write v m — Pv m and Z = PZP. 

Using Eq.(9), we transform Eq.(6) into 

Z(H + V LS )Z- X v m = E m v m , (10) 
which implies 

Z(H + V LS )Z- X = ]T £ m I w ™>< u ™ I • ( U ) 

meP 

Since E m is real (it is an eigenvalue of Eq.(5)) and 
the vectors v m are orthonormal to each other, Z(Ho + 
Vls)Z~ 1 must be Hermitian. The original problem is 
now reduced to a Hermitian model-space eigenvalue prob- 
lem 

P(H + V herm )Pv rn = E rn v m (12) 

with the Hermitian effective interaction 

V herm = Z(H + V LS )Z- X - PH P, (13) 

or equivalently 

V herm = Ze-"(H + V)e^Z- 1 - PH Q P. (14) 

To calculate Vherm, we must first have the Z transfor- 
mation. Since there are certainly many ways to construct 
Z, this generates a family of Hermitian effective interac- 
tions, all originating from Vls- For example, we can 
construct Z using the familiar Schmidt orthogonalization 
procedure, namely: 

V\ = Zuxi 

V2 = Z 2 lXl + Z22X2 

v 3 = Z 3 iXl + Z32X2 + Z 33 X3 

v± = , (15) 



with the matrix elements Z^ determined from Eq.(9). 
We denote the Hermitian effective interaction using this 
Z transformation as V sc hm- Clearly there are more than 
one such Schmidt procedures. For instance, we can use 
V2 as the starting point, which gives V2 = Z22X2, V3 = 
Z31X1 + Z32X2, an d so forth. This freedom in how the 
orthogonalization is actually achieved, gives us infinitely 
many ways to generate a Hermitian interaction, and this 
is our family of Hermitian interactions produced from 
V LS - 

We now show how some well-known Hermitization 
transformations relate to (and in fact, are special cases 
of) ours. We first look at the Okubo transformation [26| . 
From Eq.(8) we have 

(Xm I (1 +LO + U>) I Xm') = 5mm'- (16) 

It follows that an analytic choice for the Z transformation 
is 

Z = P{l+u+u) 1 l 2 P. (17) 

This leads to the Hermitian effective interaction 

V okbl = P{1+lo + uj) 1 ' 2 P{H q + Vls)P{1+oj + u J )- 1 ' 2 P 
-PH q P. (18) 

From Eqs. (8), (9), (16) and (17), it is easily seen that 
the above is equal to the Okubo Hermitian effective in- 
teraction 

V okb = P{1+L0 + LU)- 1 ' 2 {1+L0+){H Q +V) 

X(1+LU)(1+LU+LU)- 1/2 P - PH Q P (19) 

giving us an alternate expression, Eq.(18), for the Okubo 
interaction. 

There is another interesting choice for the transforma- 
tion Z. As pointed out by Andreozzi [2^, the positive 
definite operator P(l +ui + ui)P can be decomposed into 
two Cholesky matrices, namely 

P(l + uj + uj)P = PLL T P, (20) 

where L is a lower triangle Cholesky matrix, L T being its 
transpose. Since L is real and it is within the P-space, 
we have from Eq.(16) that 

Z = L T (21) 

and the corresponding Hermitian effective interaction 
from Eq.(13) is 

V cho = PL T P(H + V LS )P(L- 1 ) T P - PHqP (22) 
This is the Hermitian effective interaction of Andreozzi 
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FIG. 1: Eigenvalues of our model 5x5 interaction as a func- 
tion of the strength parameter, x. 



TABLE I: Comparison of 3x3 Hermitian (okb, cho, schm) and 
non-Hermitian (LS) effective interactions. The matrix model 
of Eq.(24) is employed with interaction strength x=0.1. 



-0.0368 0.3657 0.3976 
0.3793 0.8811 0.4482 
0.4722 0.4728 1.3999 


Vls 


-0.0358 0.3732 0.4368 
0.3732 0.8814 0.4609 
0.4368 0.4609 1.3986 


Vokb 


-0.0041 0.4063 0.5215 
0.4063 0.9020 0.5143 
0.5215 0.5143 1.3461 


Vcho 


-0.0251 0.3823 0.4698 
0.3823 0.8846 0.4695 
0.4698 0.4695 1.3846 


Vschm 



The final Hermitian effective interaction we consider 
is that of Suzuki and Okamota [27], |3J] , which is of the 
form 

V suzu = Pe- G (H + V)e G P - PH Q P (23) 

with G = tanh^ 1 ^ — oj^) and G< = — G. It has been 
shown that this interaction is the same as the Okubo 
interaction |27|. In terms of the Z transformation, it is 
readily seen that the operator e _G in Eq. (23) is equal to 
Ze~ u with Z given by Eq. (17). Thus, three well-known 
and particularly useful Hermitian effective interactions 
indeed belong to our family. 



III. MODEL CALCULATIONS 

The results from the previous section show that a fam- 
ily of Hermitian effective interactions can be derived from 
a Schmidt-type transformation of the non-Hermitian in- 
teraction Vls- We now check to see if these interac- 
tions do reproduce some of the eigenvalues of the original 
Hamiltonian, and how the effective interactions given by 
the various methods differ. In this section, we shall use 
a solvable matrix model to study these questions. Hoff- 
mann et al. |29l ] have employed a matrix model to study 
the influence of intruder states on effective interactions. 
Since we are also interested to see how intruder states 
might effect our Hermitian potentials, we use a matrix 
model of this type to study V schm , V okb and V cho , to- 
gether with their parent non-Hermitian interaction Vls- 

We employ a 5-by-5 matrix model H = Hq + xV, x 
being a strength parameter. We take Hq = {1,1,1,5,9} 



and 
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(24) 



Our P-space is chosen as the space spanned by the three 
lowest eigenstates of Hq, namely PHqP — {1,1,1}. The 
rather large diagonal matrix elements are used for V so 
that intruder states will enter as the strength param- 
eter x increases. How this happens is shown in Fig. 1, 
where the eigenvalues of H are plotted as a function of the 
strength parameter. In this figure the states are labelled 
as Pi, Pi, Pa, Qi and Qi, according to the structure of 
their wave functions. The P states are those whose wave 
functions are dominated by their P-space components, 
i.e. In contrast, the Q states are dominated by 

their Q-space components. For a weak interaction (small 
x), the lowest three states are all P states. As the in- 
teraction strength increases, states Qi and Qi decrease 
in value and states Pi and Pi increase, and we would 
expect these eigenvalues to intersect. Of course, they 
don't actually cross, but at certain interaction strength 
a Q-state "intrudes" into the P-space, becoming lower 
than the rising P-state. At x=0.5, for example, the low- 
est three states are P3, Qi and Qi, so Qi and Qi are 
intruder states in the sense that they have entered the 
P-space when the interaction is strong. We want our 
model-space effective Hamiltonian PH e ffP to reproduce 
the lowest three states of H . Thus at large x, we are re- 
quiring our PH e f fP to reproduce two intruder states. In 
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TABLE II: Comparison of 3x3 Hermitian (okb, cho, schm) and 
non-Hermitian (LS) effective interactions. The matrix model 
of Eq.(24) is employed with interaction strength x=0.55. 



-2.6246 -3.4921 -0.9479 
0.8528 1.1730 2.7893 
0.2651 0.7817 -3.0008 


V LS 


-1.7918 -0.6049 0.0615 
-0.6049 0.6264 0.7647 
0.0615 0.7647 -3.2870 


Vokb 


-1.2079 0.4599 0.6813 
0.4599 -0.4459 1.3268 
0.6813 1.3268 -2.7985 


V c ho 


-1.9123 0.3729 0.4656 
0.3729 -0.5113 1.1263 
0.4656 1.1263 -2.0288 


Vschm 



the rest of this section we study how these intruder states 
effect our potentials and what impact they have on the 
Hcrmitization procedure. 

The effective interactions are calculated using the fol- 
lowing procedures. First we calculate the wave opera- 
tor lo using the Lee-Suzuki iteration method developed 
by Andreozzi [2^|; we denote this method as the ALS 
method, and Vls is then given by Eq.(7). V sc h m is ob- 
tained from Eqs.(9), (13) and (15), and V c ho is calculated 
using Eqs.(20) and (22). For the Okubo interaction, it 
is convenient to calculate it using the method of Suzuki 

et al. [HIH Witn tllis metnod ' we nrst nnd tne 
eigenvalues and eigenfunctions defined by 



(1 + lu + lu) | a) 



a). 



(25) 



Then the Okubo Hermitian effective interaction is given 
by 



\V okb \/3)=D(a,f3) 



(26) 



with 



D(a,/3) = y/ni + 1 + J 4 + 1 



(27) 



Using these methods we can calculate Vls, V kb, Vcho, 
and Vgchm for the model potential of Eq. (24) . Two dif- 
ferent strength parameters, x — 0.1 and x = 0.55, have 
been used to see how the intruder states influence the 
effective interactions. Our results for these two param- 
eters are shown in Tables I and II. First let us inspect 
the Hermiticity of our effective interactions; clearly V Q kb, 
V c ho, and Vschm are all Hermitian, irrespective of the 



TABLE III: Comparison of eigenenergies and wavefunctions 
for the Hermitian {okb, cho, schm) and non-Hermitian (LS) 
effective interactions. The matrix model of Eq.(24) is em- 
ployed with interaction strength x=0.55, where its 5 eigenval- 
ues are -0.6714, 1.9146, 6.1702, 9.2084, and 17.4281. 





Eigenenergy 


Wavefunction 


V LS 


-0.6714 


(-0.959, 0.193, 0.207) 


Vokb 


-0.6714 


( 0.881, -0.466, -0.081) 


Vcho 


-0.6714 


(-0.528, 0.790, 0.311) 


Vschm 


-0.6714 


(-0.959, 0.193, 0.207) 




1.9146 


( 0.916, 0.397, -0.051) 




1.9146 


( 0.472, 0.853, 0.223) 




1.9146 


(-0.826, -0.563, 0.030) 




1.9146 


( 0.235, 0.951, 0.201) 




6.1702 


(-0.325, -0.179, 0.929) 




6.1702 


(-0.195, 0.932, 0.306) 




6.1702 


( 0.198, -0.241, 0.950) 




6.1702 


(-0.035, -0.234, 0.971) 



strength parameter, as they should be. The degree of 
non-Hermiticity of Vls- however, is highly dependent 
on x. In Table I, we see that Vls is only slightly non- 
Hermitian, as the largest difference in symmetric matrix 
elements is only of the order of 20%. The impact of the 
strength parameter on the non-Hermiticity of Vls can 
be seen when comparing this with table II. Here Vls is 
strongly non-Hermitian with some symmetric matrix el- 
ements differing by more than a factor of 4. Thus we see 
that when no intruder states are present (low strength 
parameter), our parent interaction is approximately Her- 
mitian, but when the intruder states enter (high strength 
parameter), our parent interaction, namely Vls, loses 
that Hermiticity in a striking manner. 

The next point to note concerning Tables I and II is the 
differences between the Hermitian effective interactions 
given by the various methods. In Table I, we see that 
the Hermitization procedures produce potentials which 
do not differ greatly from the parent potential. This, 
however, is expected since the parent potential is al- 
ready approximately Hermitian. For a high strength pa- 
rameter, we see that the resultant Hcrmitized potentials 
(V kb, Vcho, Vschm) are all indeed quite different. Thus, 
if we want Hermitian potentials which are similar, it is 
crucial that the influence of intruder states be minimal. 

Finally, we examine Table III where we show the 
eigenenergies and wavefunctions for the parent and Her- 
mitized potentials. As expected, we see that the Hermiti- 
zation procedures preserve eigenenergies. Note, however, 
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FIG. 2: The phase shifts obtained from the non-Hermitian 
L-S Viow-k are compared with those obtained from the three 
Hermitian Vi ow -k interactions in the 1 5'o partial wave channel. 



FIG. 3: The non-Hermitian Lee-Suzuki Vi ow -k interaction is 
compared with the three Hermitian Vi ow -k interactions in the 
1 Sq channel. 



the eigenfunctions of the various interactions are very dif- 
ferent, although they correspond to the same eigenvalues. 
We note also that the ground state wave function of both 
Vls and Vschm are equal to the P-space projection of the 
ground state wave function of the full-space Hamiltonian. 

IV. PHASE SHIFT EQUIVALENCE 

The non-Hermitian Vls given by the Lee-Suzuki (or 
folded diagram) method is specifically constructed to pre- 
serve the half-on-shell T-matrix T(p',p,p 2 ) [lfilllfij ; this 
interaction of course preserves the phase shift which is 
given by the fully-on-shell T-matrix T(p,p,p 2 ). It would 
be of interest to study if phase shifts are also generally 
preserved by the Hermitian interactions generated using 
the transformations described in Section II. 

Let us consider two T-matrices Ti(u>x) = V\ + 
Vig^ui^T^uji) and T 2 (lo 2 ) = V 2 + V 2 g 2 (uj 2 )T 2 (uj 2 ) , 
with the propagators <?i(cJi) = — P Hg and similarly for 
52(^2). The unperturbed state is denned by Ho, namely 
Ho\q) = q 2 \q}- The symbol P denotes the principal value 
boundary condition. These T-matrices are related by the 
well known two-potential formula 

T\ = Ti + Tl{g 2 - gx)T x + Q\{V 2 - V x )n x , (28) 

where the wave operator Q is defined by T\{uji) — 
V\£li(tjJi) and similarly for Q 2 . Applying the above rela- 
tion to the half-on-shell T-matrices in momentum space, 
we have 

(p'\Tl(p l2 )\p) = (p'lUp^p) 



+ (p'\tUp' 2 )(92(p' 2 ) ~ 3i(p 2 ))Ti(p 2 )b) 
+ (Mp')\(V2-V 1 )\Mp))- (29) 

Here the true and unperturbed wave functions are related 
by \ipi(p)) — Qi(p 2 )\p) and similarly for ip 2 . 

Using the above relation, we shall now show that the 
phase shifts of the full-space interaction V are preserved 
by the Hermitian interaction Vherm, for momentum < A. 
Let us denote the last term of Eq. (29) as D(p',p). We use 
Vherm for V 2 and V for V\ . Recall that the eigenfunction 
of (H + Vherm) is v m (see Eq.(12)) and that for H 
: II 4- V) is ^ m . We define a wave operator 

U P = |«m>(*m|. (30) 

raEP 

Then \v m ) = Up\V m ) and PV herm P = U P (H + V)U P - 
PH P. For our present case, (tp 2 (p')\ is (iy| and \ipi(p)) 
is \^f p ). Since {v v >\Up = (^f p >\, we have 

D(p',p) = (%,\(HU P -U P H)\%}, 

= (p' 2 ~p 2 )(v P >\%). (31) 

Clearly D(p,p) = 0. The second term on the right hand 
side of Eq.(29) vanishes when p 1 — p. Hence 

(p\T herm (p 2 )\ P ) = ( P \T(p 2 )\p), p < A, (32) 

where Therm is the T-matrix for (Hq + Vherm) and T for 
(Ho + V). Consequently the phase shifts of V are pre- 
served by Vherm- Recall that our T-matrices are real, 
because of the principal-value boundary conditions em- 
ployed. 
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FIG. 4: The non-Hermitian Lee-Suzuki Vi ow ~k interaction is 
compared with the three Hermitian Viow-k interactions in the 
3 Si channel. 
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fective interactions plotted with respect to momentum, using 
a model space cutoff of A = 2.0fm _1 . 



To double check this preservation, we have calcu- 
lated the phase shifts of the various Hermitian potentials 
(Vakb,V c ho and V sc hm) together with Vls- As shown in 
Fig. 2 the 1 S$ phase shifts obtained from them all agree 
with each other very well, as is the case for other par- 
tial waves. The calculations were performed with the 
CD-Bonn potential and A = 2.0/m -1 . Since the phase 
shifts of Vls are, by construction, the same as those of 
the full-space potential V, the Hermitian potentials pre- 
serve the phase shifts of V. It has been pointed out Q 
that the Okubo Hermitian potential preserves the phase 
shifts, but we have found that there is a family of Hermi- 
tian potentials, including Okubo, which all preserve the 
phase shifts up to the decimation scale A. 



V. HERMITIAN LOW-MOMENTUM 
INTERACTIONS 

A main purpose of having a low- momentum nucleon in- 
teraction Vi ow -k is to use it in nuclear many body prob- 
lems such as the shell model nuclear structure calcula- 
tions. As we have seen, however, there is a family of 
phase-shift equivalent Hermitian interactions. How dif- 
ferent are they? Which one should one use for nuclear 
structure calculations? Will these Hermitian effective in- 
teractions give rise to the same physical properties? It is 
these questions we now seek to answer. 

We have calculated the Hermitian Vi ow -kS corre- 
sponding to Vokb, V c ho and V schm using several NN po- 
tentials. Calculations for the non-Hermitian V/ oto _fc cor- 
responding to Vls were also performed. In Fig. 3 we 



compare the results for the 1 Sq channel, obtained with 
the CD-Bonn potential and A = 2.0 fm~ l . Clearly the 
Hermitian interactions are all quite similar to each other 
and to the parent non-Hermitian potential Vls- In Fig- 
4, we show a similar plot for the 3 Si channel. Again, with 
the exception of very low momentum, the Hermitian po- 
tentials are all nearly identical to the parent Vls- It is of 
interest that at very low momentum, the V sc h m matrix 
elements are slightly more attractive than the others. We 
note that the Hermitian effective interactions all preserve 
the deuteron binding energy (2.225 MeV). In addition, 
they are all phase shift equivalent up to the decimation 
scale A, as illustrated in Fig. 2. 

Despite the similarities, there is, however, an addi- 
tional degree of preservation that V sc h m satisfies but the 
other two Hermitian interactions don't. By construction, 
the deuteron ground-state wave function given by V sc hm 
is exactly equal to the P-space projection of the wave 
function of V, which is not true for V kb and V c ho- This 
additional preservation is worth studying further. We 
refer to Figures 5 and 6, where we plot the S and D- 
state deuteron wavefunctions for a cutoff of A = 2.0fm -1 . 
Overall, the agreement is very good between all poten- 
tials, not just Vls and V sc h m , which is to be expected 
considering that the interactions themselves are approx- 
imately the same. The fact that V sc h m gives exactly the 
same wavefunction as Vls can be seen from the D-state 
probability of the deuteron for each interaction, which we 
list for convenience in Fig. 6. Whereas the Pd's for the 
effective interactions are close, they are exactly equal for 
Vls and V sc h m - This exact preservation of the ground 
state wavefunction presents us with an extra constraint 
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TABLE IV: A comparison of shell model relative matrix element calculated using Vls with those calculated using the Hermitian 
interactions. The oscillator length parameter b is given by b — -J -^j, in units of fm. The matrix elements are in units of 
MeV. 
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that might be useful in deciding which Hermitian poten- 
tial to use. 

To ensure our Hermitian potentials will be useful for 
nuclear structure calculations, we examine in Table IV 
some shell model matrix elements calculated with Vi ow -k 
corresponding to V ok b, V cho , V schm and V L s- As seen they 
are all virtually identical, as would be expected owing to 
the similarity between the potentials themselves. This 
good agreement is a desirable result, as it implies that 
shell model calculations will not be sensitive to which 
interaction one employs. 

We have seen that the Hermitian potentials generated 
above are all approximately the same, and we would like 
to offer an explanation as to why this is so. Although 
the Vi ow -k corresponding to Vls is non-Hermitian, we 
emphasize that it is only slightly so. In reference to our 
model study of Section 2, this corresponds to the situa- 
tion with a small strength parameter, and thus it would 
not be surprising to see the Hermitian Vi ow -k potentials 
are so similar to the parent Vi ow -k- For example, this is 
especially transparent in the case of the Okubo interac- 
tion of Eq. (18) or (19), as it can be written as [27Ll35| 



(alVok. 



kb 



H\{V ls 



V S= (a\-(Vr,s - Vl s )\p) . 



v£s)\P) 



(33) 



This tells us that V Q kb can be well approximated by Vls-, 
if Vls is only slightly non-Hermitian. In fact in this case 
V kb is very accurately reproduced by the simple average 
(V LS + vl s )/2. 

We note that the D-state probabilities given in Fig. 6 
for the various effective interactions are significantly dif- 
ferent from that of the bare potential. This is due to 
the fact that the effective interactions are subjected to a 
momentum cutoff at 2.0 fm -1 , while the bare CD Bonn 
potential extends far beyond that. This can clearly be 
seen from the figure where a large portion of the wave- 
function is simply cut off. While Pjj is not an observable, 
it has theoretical relevance-it is an important character- 
istic of modern nucleon-nucleon potentials. 

According to the MEEFT prescription, it is necessary 
to impose the cutoff at the limit of experimental data 
(in this case the limit of NN scattering experiments at 
2.0 fm -1 ), and it is at this cutoff that the effective inter- 
actions from all the bare potentials become nearly iden- 
tical. Raising the cutoff would erode this approximate 
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FIG. 6: D-state deuteron wavefunctions for the bare and ef- 
fective interactions plotted with respect to momentum, using 
a model space cutoff of A = 2.0fm~ 1 . 



FIG. 7: D-state deuteron wavefunctions for the bare and ef- 
fective interactions plotted with respect to momentum, using 
a model space cutoff of A = 3.0fm _ . 



uniqueness and model independence. But at the above 
cutoff, we can not preserve Pn ■ To illustrate, we refer to 
Figure 7 where we plot the D-state wavefunctions with a 
cutoff of A = 3.0 fm -1 . It shows that much more of the 
wavefunction is contained in the region below the cutoff, 
and as a result, the Pd's for the effective interactions are 
much closer to that of the bare potential. If we increase 
the cutoff to 4.0fm _1 , they are almost exact. How to re- 
solve this disparity in how to choose the cutoff warrants 
further study. 

VI. SUMMARY AND CONCLUSION 

We have studied a general method for deriving low- 
momentum NN interactions which are phase shift equiv- 
alent. By integrating out the high momentum compo- 
nents of a realistic NN potential such as the CD-Bonn 
potential, the Lee-Suzuki or folded-diagram method is 
employed to derive a parent low-momentum NN poten- 
tial. This potential preserves the deuteron binding en- 
ergy and phase shifts up to cutoff scale A. In addition 
it preserves the half-on-shell T-matrix up to the same 
scale. This Lee-Suzuki low-momentum NN interaction 
is not Hermitian, and further transformation is needed 
to obtain low-momentum interactions which are Hermi- 
tian. We have shown how to construct a family of such 
an interaction using the Schmidt orthogonalization pro- 
cedure, and we have seen that two existing Hermitization 
schemes, namely Okubo and Andreozzi, are in fact spe- 
cial cases of our general process. We have shown that all 
the Hermitian interactions so generated are phase shift 
equivalent, all reproducing empirical phase shifts up to 



scale A. These potentials also preserve the deuteron bind- 
ing energy. 

Through an analysis of this procedure using a solvable 
matrix model we have seen some interesting properties of 
Schmidt transformation method. In particular, with the 
entrance of intruder states, the parent potential can be- 
come highly non-Hermitian, and that the Hermitian po- 
tentials can deviate largely from each other and from the 
parent potential. It is fortunate that such deviations are 
not present for low-momentum NN potentials, for cutoff 
momentum A ~ 2.0 fm . This is mainly because our 
parent Vi ow ~k is only slightly non-Hermitian, and as a re- 
sult, the Hermitian low-momentum nucleon interactions 
generated from our orthogonalization procedure are all 
close to each other and close to the parent Vi ow -k- Shell 
model matrix elements of the low-momentum nucleon in- 
teractions are found to be approximately independent of 
the starting Hermitian potential, indicating the useful- 
ness of our procedure in nuclear many body calculations. 
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